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Abstract. We present a proposal and a feasibility study for the creation and 
quantum state tomography of a single polariton state of an atomic ensemble. 
The collective non-classical and non-Gaussian state of the ensemble is generated 
by detection of a single forward scattered photon. The state is subsequently 
characterized by atomic state tomography performed using strong dispersive light- 
atoms interaction followed by a homodyne measurement on the transmitted light. 
The proposal is backed by preliminary experimental results showing projection 
noise limited sensitivity and a simulation demonstrating the feasibility of the 
proposed method for detection of a non-classical and non-Gaussian state of the 
mesoscopic atomic ensemble. This work represents the first attempt of hybrid 
discrete-continuous variable quantum state processing with atomic ensembles. 
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1. Introduction 

Atomic ensembles have emerged as a highly efficient medium for light-matter quantum 
interfaces [1]. Mapping of squeezed states of light onto an ensemble [2, 3], generation 
and retrieval of single excitations [4, 5, 6], entanglement of two ensembles [7, 8, 9] 
and quantum sensing [10, 11, 12] have all been demonstrated experimentally. The 
development of interfaces between light and atomic ensembles have so far followed 
two separate paths which until now have been pursued independently. Operating with 
either discrete excitations and single photon detection or with continuous Gaussian 
states and homodyne measurements. Within the latter approach tomography of spin 
squeezed atomic states by quantum non-demolition (QND) interactions with light has 
been developed [13, 14, 15]. 

In close analogy to methods used in photonic systems [16, 17], the hybrid approach 
of combining discrete excitations, such as number states, and measurements in the 
continuous variable domain opens up new venues in quantum state engineering with 
atomic ensembles. Among them are hybrid quantum repeaters and generation of 
Schrddinger cat states [18], weak quantum measurements [19] and heralded quantum 
amplification for precision measurements [20]. 

In the following, we report on the first steps towards bridging the gap between 
the discrete and continuous variable approaches with atomic ensembles. Starting 
from a collective photon scattering as in Duan et al. [21], we present a method to 
produce and directly characterize a highly non-classical (negative Wigner function) and 
non-Gaussian collective (entangled) atomic state. Using interference between different 
atomic spin-wave modes allows us to demonstrate its properties. In particular, we will 
demonstrate how a single atomic excitation changes the projection noise statistics of a 
macroscopic atomic cloud containing over 100 000 atoms fundamentally. 

A central point in our proposal is that both the creation and characterization 
is performed directly in the atomic ensemble. This has several advantages compared 
to other methods: First, the direct creation inside an atomic ensemble with a long 
coherence time results in an heralded state that is readily available for on-demand use. 
Secondly, in comparison to protocols where the quantum state of the atomic ensemble is 
created by mapping of a photonic state on to the ensemble [5, 22, 23] our method does 
not suffer from the inherent loss mechanism in the transfer of a quantum state between 
light and atoms. Third in in comparison to converting the single atomic excitation 
back into into the optical domain for further characterization [24], preserving the 
state inside the memory presents perspectives for further manipulation, e.g. improved 
heralding scaling by shelving the excitation into a third atomic level, Gaussian state 
manipulation by dispersive QND-measurements and feedback, as well as linear optics 
equivalent operations between different spin-wave modes. 

2. Generation of a single polariton state 

We consider a system of N a spin-1/2 particles described by the states |f) and \l) which 
in the experiment are the two hyperfine ground states of a Cesium (Cs) atom, also 
known as the clock-levels. Initially all atoms are prepared in the |f) state, see figure la, 
such that the quantum state of the ensemble can be written as the product state 

N a 

l*o) = ®|t), = |tt..-tt)- (1) 

Z=l 
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Figure 1. The atomic levels and the collective Bloch sphere (as described by 
Dowling et al.[25]) at different stages of the protocol (a) Atoms are prepared in 
l^o) and a weak blue detuned excitation pulse is applied (b) detection of a forward 
scattered Raman photon signals the generation of a single collective excitation in 
the ensemble (c) a microwave tt/2 pulse rotates the state l^i) into the equatorial 
plane of the Bloch sphere, the insert shows a plot of the non-Gaussian (Gaussian) 
marginal distribution (probability density of J z -measurement outcomes) of the 
l^f'l) (|^o)) state after a rotation to the equatorial plane in blue (orange). 



Next, a weak off- resonant excitation pulse is sent into the ensemble, and conditioned 
on the detection of one photon scattered forward via a Raman process, one spin flip 
occurs in the collective polariton state of the ensemble. This detection event heralds 
the preparation of a single excitation of a collective atomic zero-transverse-momentum 
spin wave - a single polariton, see figure lb. The success probability of this forward 
scattering is kept low, such that the probability to forward-scatter two or more photons 
is negligible. The state of the ensemble [21] becomes 

i N ° 

I*!) = at |* ) = ]T |tt . . . t I t-.-tt), (2) 



1 1 l-th atom 
where we have defined the symmetric creation operator 

N. 



at 



= X>Mtlr ( 3 ) 



=1 

To describe the atomic ensemble we define pseudo spin operators for each individual 



atom 



j^i(|l><t| + |t)Ul), (4) 

^=-i(i4)(ti-it)ai), (5) 
h = l(n)tt\-\i){i\), (6) 
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which obey angular-momentum like commutation relations. From these we define 
symmetric (under particle exchange) collective ensemble operators by summing the 
individual atomic operators as 

Ji=^2ii ] for i = x,y,z. (7) 
1=1 

where N a is the total number of atoms in the ensemble and jf^ acts on the Z-th atom. 

A subsequent microwave 7r/2-pulse rotates this collective entangled state, see 
figure lc, such that it becomes an even superposition of |f) and \\.) given as 



|*i> =Ar/al*i) (8) 
1 ^ 

^ •+...+ - + ... + +), (9) 



l-th atom 



where the microwave it/2 pulse is described by the operator 



N a 



i2,r/2 = 2^(l+)l<t| l + |-> I UI I ), (10) 
1=1 

and the coherent superposition states are 

w-K^. ("I 

To show that the state has a non-Gaussian population-difference between the two 
hyperfine levels we calculate the probability to find n atoms in the ||) state and the 
remaining (iV a — n) atoms in the |f) state. This is done by calculating the overlap 
between the single excitation state, |*i), and the general n-th excited state given by 

/AT \-l/2 



permutations 



find 



N a — >00 n / Z „2 

« 2,y_ a-, (14) 
whereas the standard coherent spin state (CSS) \^' ) has 

K*nK>| 2 =2-»-( N °) (15) 



irN n . 



(16) 



with cc = y \ n 2)' From this it is clear that, conditioned on whether 

a single photon is detected or not, the probability distributions for outcomes of 
J y - and J z -measurements will be profoundly different, see insert on figure lc, i.e. 
the respective probability for a J z -measurement to give (n — N a /2) as outcome is 
P(J z =n-N a /2) = \{* n \% orl )| 2 . 
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In the ideal case of perfect tomography of a pure single polariton excitation the 
population difference follows the marginal distribution of a Wigner function with a 
single excitation, a n = 1 Fock-state, given by equation (14). In a realistic case of finite 
efficiency of both the state preparation and state detection we show that the marginal 
distribution still retains its non-classical and non-Gaussian features, as discussed in 
detail in section 3. 

Before we continue to the experimental implementation we introduce the effective 
quantum efficiency of the quadrature measurement, e. This is equivalent to the 
efficiency of the state tomography which is based on the quadrature measurements. 
These are done with a dispersive QND-probing method introduced in [26] and discussed 
in detail in [14, 27]. We use two probes, each of which interacts with the |f) (\D) state 
exclusively and which therefore does not cause Raman transitions to the other clock 
state \l) (|t)), respectively. Each of these probes measures the atomic state dependent 
optical phase shift, which we then convert into quadrature values [28]. By a noise 
scaling analysis we discriminate light shot noise and technical noise against the atomic 
noise. With this we quantify the efficiency as 

var (CSS(JV a )) - var (CSS(A a = 0)) 
E var(CSS(iV a )) ' ( ' 

Here CSS(iV a ) denotes the differential phase shift measured on a coherent spin state 
l^o) with N a atoms. Note that var(CSS(iV a = 0)) (the optical phase fluctuations 
measured in the absence of atoms) contains both shot noise and technical noise. The 
contribution of atomic noise to var (CSS(iV a )) scales quadratically with the number 
of photons in the probe, whereas the shot noise contribution only scales linearly [1]. 
For tomography applications it is not necessary to preserve the non-destructiveness of 
the QND probing - that is to preserve the coherence between the states |f) and |T). 
Therefore, neglecting technical noise, a probe of unlimited strength can in principle be 
used, and the quantum efficiency of the tomographic measurement will approach unity. 

3. Implementation 

As described in section 2, the Wigner function of the single polariton state has both 
a non-Gaussian marginal distribution and a negativity and is therefore highly non- 
classical. Over the past years several experiments have shown atomic state tomography 
[10, 13, 27, 29] in which the characterization of such a state could be performed. In 
the following we will present a detailed discussion on how a single polariton state can 
be created and characterized. We will focus on the specific experimental configuration 
described in [14, 27]. 

3.1. Experimental considerations 

The experiment is based on an ensemble of approximately 10 5 Cs atoms. We consider 
the two- level system formed by the clock levels, i.e. \\) = \F = 3, mp = 0) and 
|t) = \F = A, nip = 0). As the excited state, used to couple the two clock states, 
we use the \e) = \F' = 4, m! F = +1) state of the 852 nm D2 transition. The atomic 
ensemble is situated in a far-off resonant dipole trap which is placed in one arm of a 
Mach-Zehnder interferometer (MZI), see figure 2. To obtain the probability distribution 
of J z -measurement outcomes, we detect the atomic state dependent optical phase-shift 
imprinted on light passing through the atomic ensemble using a dual-color dispersive 
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Figure 2. Sketch of the experimental setup: The atomic ensemble is trapped in 
a far off-resonant dipole trap (not shown). A single polariton state is generated 
by sending a weak excitation pulse through the atomic cloud. The detection 
of a forward scattered single photon heralds the creation of the desired state. 
Polarization and frequency filtering is applied to avoid the detection of undesired 
excitation photons using a polarizing beam splitter and two Fabry-Perot filter 
cavities. To characterize the atomic state we measure the atomic state dependent 
optical phase shift imprinted on a light pulse by the atoms. 



QND measurement as described in section 2. To create the single polariton state given 
by equation (2) we propose the following procedure; initially all atoms are prepared in 
the |t) state by a combination of microwave pulses, optical pumping and purification 
pulses [30]. The collective excitation is created by a weak excitation pulse, blue detuned 
by A ss |r from the |f) — > \e) transition. The detuning is chosen such that the optical 
depth is less than unity, (assumed an on resonant optical depth of OD « 10) which 
means that all atoms will have about the same chance to scatter a photon. If this 
condition is fulfilled the excitation will be collective, i.e. shared almost uniformly by 
the whole ensemble. 

To obtain the best possible spatial overlap between the probe beams and the 
excitation beam with the atomic ensemble all beams arrive at the MZI input via 
the same optical fiber, whose output is carefully overlapped with the atomic sample. 
This nearly perfect overlap of the light beams obtained by using the same optical 
fiber comes with the constraint that the polarization of the excitation beam and the 
probe beams are now identical. With the quantization axis (in z-direction, defined 
by the magnetic bias-field) chosen orthogonal to the propagation direction of the 
beams (y-direction) , we can only address the atoms with it- or x-polarized light 
(i.e. linearly polarized orthogonal to 7r-polarization). In order to obtain the highest 
efficiency of the atomic tomography it is preferable to use 7r-polarized light [14]. This 
combined with the complicated multi-level structure of the Cs atoms puts several 
constraints on the excitation and probe geometry and necessitates the application of 
a strong magnetic bias field as discussed in the following. Because of the symmetry 
of the fiif = atomic wave functions, anti-Stokes photons on the Raman transition 
|t) — > \F' = 4) — > \l) cannot be emitted into the forward direction in the absence 
of an external magnetic field. Since for both x- and 7r-polarized excitation beams, 
the Clebsch-Gordan coefficients corresponding to the relevant dipole matrix elements 
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Figure 3. Cs level structure and excitation scheme: upwards arrows indicate 
the excitation light and the downward arrows indicate decay channels. Light 
corresponding to the dashed arrows is filtered out by either polarization (blue) 
or frequency (orange). Solid arrows denote decays that cannot be filtered out. 
Note that the decay to \F = 3, rap = 2) leads to the unavoidable inefficiency of 
the state generation as discussed in the text. 

are such, that excitation can only occur into superpositions of excited state Zeeman- 
sublevels for which the same-polarization decay into the respective other \F, xtif = 0) 
hyperfine ground state interferes destructively in the absence of magnetic fields. The 
applied bias field is approximately 20 Gauss which shifts the \F' = &,m'p = —1) state 
out of resonance by several line widths, see figure 3. This choice combined with the fact 
that the |f) — > \F' = 4, m' F = 0) transition is forbidden by selection rules allows us to 
achieve the required selective excitation to the |e) state. From the |e) state the atom 
can undergo spontaneous emission through six possible decay channels, see figure 3. 

We now apply different filtering methods such that we only detect photons 
corresponding to an atom decaying via the |e) — s- \\) transition, since this projects the 
ensemble into the the single polariton state given by equation (2). Photons with tt- 
polarization are suppressed by a polarizing beam splitter. The harder task is to filter out 
photons that are x-polarized and have a frequency corresponding to the |e) — > \F = 4) 
manifold. Photons with these properties originate from both the excitation beam and 
the spontaneous decay. We thus need to reject photons with a frequency difference 
of approximately 9 GHz. This is done with two cascaded Fabry-Perot filter cavities, 
see figure 2. This way we filter out all photons except the ones originating from 
the undesired decay channel |e) — » \F = 3, mp = 2). This decay channel produces 
forward scattered photons which have the same polarization as the desired photon, 
the frequency difference between them is in the MHz range (the Zeeman splitting of 
the ground state), which makes it experimentally very hard to selectively reject them. 
Note that the branching ratios favor the preferred decay by a factor of « 4. 

With the single photon detected, we proceed in a fashion similar to J. Appel et 
al. [27]. We apply a 7r/2-pulse and perform the dispersive QND probing of the atomic 
population difference in ||) and ||). Repeating this procedure several thousand times 
we characterize the distribution of J z -measurements for the prepared atomic state. 
Note that due to the rotational symmetry of the Wigner function of both the single 
polariton- and the coherent spin-state the ir/2 rotation can be performed without 
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Figure 4. The measured variance of the probe phase as a function of atom number, 
with different noise contributions discriminated by noise scaling analysis. Blue 
line: general quadratic fit to data (blue stars with statistical errors). The atomic 
projection noise (green area) scales linearly with N a . Technical fluctuations (red 
areas) originate both from noise in the QND measurement as well as imperfections 
in the CSS preparation. 

paying attention to the phase of the microwave pulse (which defines the axis of rotation 
on the Bloch sphere). 

3.2. Test of tomography noise performance 

The major experimental challenge is concerning the characterization of the atomic 
state via quantum tomography, since this requires the measurement of the population 
difference (J z ) in the experimental apparatus to be projection noise limited. We have 
checked our noise sensitivity in the experimental configuration with a bias-field of 
B ss 20 Gauss, corresponding to a splitting between the |e) and \F' = 4, m' F = —1) 
states of several line widths. 

We start by preparing a CSS of the atomic ensemble (a product state of each 
atom being in |+)) [27] and then measure the population difference between the two 
clock-levels |f) and Repeating this several thousands times and performing a 
noise scaling analysis using the same method as in [27], we attribute the measured 
fluctuations to different origins, see figure 4. We see a predominantly linear dependence 
of the atomic noise on the atom number (green area) - a clear signature of the required 
projection noise limited sensitivity. This certifies sufficient performance of the initial 
state preparation, the quality of our microwave source and the long term stability of 
our setup. 

The next part of the experiment is to detect the single photon, distinguishing and 
filtering out photons originating from the excitation beam and undesired decay channels. 
Detection of any of these photons would result in false positive "clicks" , i.e. we detect 
a photon but the atomic ensemble is not prepared in the desired single excitation 
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Table 1. Probability for photo-counts of different origins. We have only kept 
the most dominant processes, i.e. probabilities greater then 0.05%. 



Origin of photo-count 


Created state 


Probability (%) 


Dark counts 


l*0> 


8.7 


Leakage of excitation pulse 


1*0} 


0.6 


Decay via |e) — > \F = 3, rap = 2} 


l*0> 


19.1 


Decay via |e) — > |4) 


l*l> 


71.5 



state | \& i ) . Such "clicks" could be due to dark counts, leakage of the probe-, trap- 
or excitation-beams. Since the state generation is based on the detection of a single 
photon with a low success probability, false positives essentially mix in realizations 
prepared in the vacuum state (zero polaritons). 

To describe this, we consider the actual state prepared (after the tt/2 microwave 
pulse) conditioned on a "click" as a classical mixture of the single excitation state 
which is obtained with a probability p, and the CSS, \^q), obtained with a probability 
(1-p): 

p = p\%)(*[\ + (l-p)\%)(%\ (18) 

In order to show that the outlined procedure is suitable for our experimental apparatus 
we perform an analysis based on a simulation using the relevant experimental 
parameters. 

4. Analysis 

To quantify the non-Gaussian character of the prepared state we calculate the expected 
probability distributions for our experimental parameters and perform a simulation 
similar to Dubost et al. [31]. 

In order to estimate the purity of the single polariton excitation (the probability 
p introduced in the previous section), we assume that the excitation pulse is so weak 
that the probability of producing a forward-scattered photon is 5 %. To reach this 
scattering probability we require 1.35 • 10 4 photons in a 10 fis long excitation pulse. 
The polarization filtering is done with a polarizing beam splitter cube with a rejection 
of 1 : 7 • 10 3 . For the frequency filtering we obtain rejections of 1 : 5 • 10 7 by using two 
cascaded filter cavities, each with a transmission of 80%. The mode-matching overlap 
of the single photon detection mode and the atomic state tomography mode is taken 
to be 75%. With this we estimate the probabilities of the photo-count originating 
from a specific decay channel, a dark count of the detector or from the leakage of other 
light sources (e.g. the excitation pulse) in the experiment. The probabilities of the 
dominating processes (those exceeding 0.05%) are shown in table 1. 

From the probability given in table 1 we calculate the probability distributions for 
^-measurements. The created state is modeled as a statistical mixture of a CSS and a 
ideal, noise free single excitation state weighted with their corresponding probabilities, 
see equation (18), a plot is shown in the insert on figure 5. The efficiency of the 
atomic state tomography is reduced due to the extra noise added by the readout 
procedure: We use light as the meter-system and therefore the shot noise of light 
will also have a contribution. Following the approach outlined in [32] we model this 
independent Gaussian noise as an extra admixture of the vacuum state to the state to 
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Figure 5. Calculated distributions for j z -measurements outcomes of the created 
single polariton state (blue) and the vacuum state (orange), taking into account 
the finite tomography efficiency. The histogram is based on 1000 randomly drawn 
samples from each distribution. Insert: Marginal distributions not including the 
effects of a non-unity detection efficiency. Note that the less-than-unity tomography 
efficiency washes away the non-Gaussian feature, making it hard to distinguish 
the two histograms by eye. 

be analyzed. In effect this reduces the non-Gaussianity of the J z distribution detected 
by the dispersive measurements as shown in the main panel of figure 5. 

As pointed out in [31], the detection of a non-Gaussian marginal distribution 
is an experimental challenge. One of the main reasons for this is that we can only 
estimate the underlying probability distributions by a finite number of experimentally 
acquired samples. To illustrate this, we draw 1000 samples from the distributions for 
the dispersive J z -measurements. These samples are binned and plotted as a histogram 
in figure 5. It is clear from the histograms that distinguishing the two sets of sampled 
data and thus detecting the non-Gaussianity and possibly the non-classicality of the 
quantum state is experimentally challenging. Several methods have been developed in 
order to quantify the non-classicality [33, 34], we note that these methods define the 
term "non-classical" in different ways. The criterion by Kot et al. does not implicitly 
assume quantum mechanics, whereas the less stringent Vogel criterion does: A state 
described by equation (18) is non-classical for all p according to the Vogel criterion 
[35], whereas only states with p > 0.5 are non-classical by the definition of Kot. et al. 

In the following we present a simulation based on the ideas in [31] where the 
non-Gaussianity is quantified in terms of the statistical cumulants. From the known 
probability distributions we draw random samples and use these to calculate the second 
and fourth cumulants, K2 = M2 an d K4 = (14 — 3/x|, for varying sample sizes, Hi denotes 
the i-th statistical moment. Repeating this 1000 times allows us to calculate the 
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Figure 6. Plot of the cumulants K2 and function of the sample size. 

Dotted points correspond to data from one realization simulation. 



corresponding standard deviations, see figure 6. We start by considering the second 
order cumulant, which is just the variance of the distribution: From figure 6a we note 
that only around 250 samples are needed to distinguish between the distribution of the 
single polariton state and the vacuum (CSS). We thus expect that careful investigation 
of K2 will allow us to clearly classify our data as originating from cither the expected 
single polariton state or the CSS under the assumptions that these are the only possible 
distributions. In order to verify that our state is non-Gaussian we need to observe 
a non-zero forth cumulant K4. We expect to require approximately 500 samples to 
clearly distinguish the observed probability distribution form the Gaussian one of 
the CSS with K4 = 0. The performed simulations are based on realistic parameters, 
many already experimentally verified, and we can thus conclude that our proposed 
implementation is experimentally feasible. 



5. Conclusion 



In this paper we have introduced a hybrid method based on discrete excitations and 
continuous measurements which allows us to generate and characterize a single polariton 
state of an atomic ensemble. This state has non-Gaussian marginal distributions of 
the Wigner function and is non-classical. We have presented a detailed proposal for 
the experimental creation and detection of the single polariton state. The proposal is 
backed by a simulation for experimental valid parameters together with preliminary 
results showing the feasibility of the proposal. This work is a step towards implementing 
a hybrid approach to quantum information processing with atomic memories as well 
as a contribution to the ongoing research into the foundations of quantum mechanics. 
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